function plot_h_refine_sample
global V T E ET TE posT posV posE;

%1. refine case 1:
V = [0,0;1,0;1/2,1/2;1/2,-1/2];
T = [1,2,3;4,2,1];
E = [1,3;2,3;1,2;1,4;2,1];
posT = 2; posV=4;posE = 5;
TE=[1,2,3;3,4,5];
ET=[1,0;1,0;1,2;2,0;2,0];
% refine
[v_eg,eg,eg1,eg2] = refine_me(1);
refine_nei_green(2,v_eg,eg1,eg2,eg); 
%plot
figure; plot_t(V,T);
Bdr=[1;2;3;1;4;2];
hold on; plot(V(Bdr,1),V(Bdr,2),'r');plot(V(:,1),V(:,2),'o');hold off;

%1. refine case 2:
V = [0,0;1,0;1/2,1/2;0.9,-2];
T = [1,2,3;4,2,1];
E = [1,3;2,3;1,2;1,4;2,1];
posT = 2; posV=4;posE = 5;
TE=[1,2,3;3,4,5];
ET=[1,0;1,0;1,2;2,0;2,0];
% refine
[v_eg,eg,eg1,eg2] = refine_me(1);
refine_nei_green(2,v_eg,eg1,eg2,eg); 
%plot
figure; plot_t(V,T);
Bdr=[1;2;3;1;4;2];
hold on; plot(V(Bdr,1),V(Bdr,2),'r'); plot(V(:,1),V(:,2),'o');hold off;

%1. refine case 3:
V = [0,0;1,0;1/2,1/2;0.1,-2];
T = [1,2,3;4,2,1];
E = [1,3;2,3;1,2;1,4;2,1];
posT = 2; posV=4;posE = 5;
TE=[1,2,3;3,4,5];
ET=[1,0;1,0;1,2;2,0;2,0];
% refine
[v_eg,eg,eg1,eg2] = refine_me(1);
refine_nei_green(2,v_eg,eg1,eg2,eg); 
%plot
figure; plot_t(V,T);
Bdr=[1;2;3;1;4;2];
hold on; plot(V(Bdr,1),V(Bdr,2),'r');plot(V(:,1),V(:,2),'o');hold off;